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Abstract 

Motivated by the recent article of P. Shea et al. [Am. J. Phys. 77 (6), 2009] we examine the 
exactly solvable problem of two harmonically trapped ultra-cold bosonic atoms interacting via a 
short range potential in one and two dimensions. A straightforward application in one dimension 
shows that the energy spectrum is universal, provided that the range of the potential is much 
smaller than the oscillator length, in addition to clearly illustrating why regularization is not 
required in the limit of zero range. The two dimensional problem is less trivial, requiring a more 
careful treatment as compared to the one dimensional case. Our two dimensional analysis likewise 
reveals that the low-energy physics is also universal, in addition to providing a simple method for 
obtaining the appropriately regularized two dimensional pseudopotential. 
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I. INTRODUCTION 

Recently, P. Shea et al. [1] have discussed the problem of two bosonic atoms interacting 
via a short range potential and trapped in a three dimensional (3D) spherically symmet- 
ric harmonic oscillator potential. Their work showed that the low energy properties of the 
3D system are universal, irrespective of the shape of the potential, provided the range is 
much smaller than the oscillator length. In addition, they developed the appropriate pseu- 
dopotential for a zero-range interaction in a manner accessible to undergraduate students 
of physics with only an elementary knowledge of quantum mechanics and scattering theory. 
Specifically, no prior knowledge of self-adjoint extensions, renormalization techniques, or 
dimensional regularization schemes are required to obtain equivalent results. 

In this article, we present details of the application of the techniques presented in Ref. [l| 
to both one and two dimensional systems. While the one dimensional (ID) problem proves 
to be quite accessible, the two dimensional (2D) analogue turns out to be rather subtle. 
The current interest in the physics of low-dimensional cold atom systems [2j provides addi- 
tional impetus for the results presented in this paper. In particular, it is now conceivable 
that analogous experiments to those performed in Ref. Q can also be carried out on low 
dimensional systems, in which case the universal aspects of the spectra derived here may be 
experimentally verified. 

The plan for our paper is as follows. In section II, we derive the pseudopotentials appro- 
priate for a two-body interaction in the limit of zero range for both the ID and 2D systems. 
Our approach clearly illustrates the concept of the pseudopotential in the form of a regu- 
larized Dirac delta function while avoiding technical discussions about the self-adjointness 
of the two-body Hamiltonan, or of regularization operators required to ensure the Hamil- 
tonian's self-adjoint property. Our result for the ID pseudopotential scenario agrees with 
the literature, whereas our 2D pseudopotential is ostensibly different from earlier published 
results. Nevertheless, we argue that our 2D pseudopotential is operationally equivalent pro- 
vided that it acts upon the appropriate two body wave function. In section III, we show 
that the energy spectra in both the ID and 2D systems are universal and independent of 
the details of the pseudopotential provided that the range of the interaction is much smaller 
than the oscillator length. The universal properties for the energy spectrum we find in 
both ID and 2D are not well-known in the literature. [4J Our analysis also provides a sharp 
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contrast to the results presented in Ref. (5( where the explicit properties of the regularized 
Dirac delta function are needed in order to obtain the spectrum. In section IV, we finish 
with some concluding remarks and suggestions for future research in this area. 



II. THE ZERO-RANGE PSEUDOPOTENTIAL 



One dimensional treatment 



We first consider a free system of two identical bosons, each of mass M, interacting via a 
short-range symmetric potential in one dimension. In the relative coordinate, r = |r*i — r^l > 
0, the s-state asymptotic scattering wave function is given by 

ij){x) ~ cos(£x + 5{k)) (r > b), (1) 

where b is the range of the interaction potential and r = \x\. Barlette et al. 6| have already 
provided us with the s-wave effective range expansion, relating the phase shift, 5(h), to the 
effective range, r , and scattering length, a, viz., 

k ttm(5(k)) = - + K k 2 + 0(k 4 ) , (2) 

where the higher order terms are shape dependent. [l| The effective range, r , is related to 
the range, b, in such a way that as b — > 0, r n — > 0. Thus, in the limit of zero-range, Eq. (J2]) 
reduces to fctan(5(/c)) = -. Following Ref. [l|, we extrapolate these results to bound states 
for positive a, where the S'-matrix has poles at cot S = i. Utilizing the zero-range limit of 
Eq. (j2J), we immediately obtain k = which gives the bound state wave function, 

^{x) = e'-- (r > 0), (3) 

with binding energy E = h 2 k 2 /M = —jf^- Equation (3) is exact for a zero-range potential 
and holds for b ^ provided the size of the bound two-body system is much larger than 
the range of the potential responsible for the binding. Under these conditions, it is possible 
to construct an effective potential that reproduces the shape-independent results we have 
just obtained. To this end, we consider the ID Laplacian of Eq. (3) with respect to the 
argument, and recall that r = \x\. Performing this operation gives 

d 2 ?p(x) 1 _m f d\x\\ 2 1 _nd 2 \x 



<li 2 a 2 \ dx ) a 6 dx 2 ^ ^ 



Using the relations (^fj = 1 and = 25 ld (x) in Eq. (jlj) yields 



h 2 d 2 ij(x) 2h\ u . _N ft 2 _N _ 

-M^-WJ {x)e " =-M^ 6 (5) 

which, upon the substitutions E = — and ip(x) = e~", reduces to 



H 2 d 2 tp{x) 2ti 
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■5 ld (x)ip(x) = EiP(x). (6) 



M dx 2 Ma 

Equation ([6]) it nothing more than the time-independent Schrodinger equation from which 



we see that the pseudopotential 



V p \ D (x) = ~^(x), (7) 



reproduces the earlier results for the wave function and binding energy for any zero-range 



potential with a given scattering 
agreement with the literature, p, LZ 



ength. This zero-range pseudopotential is in complete 
, Is] and naturally illustrates why no regularization of the 



Dirac delta function is required in one dimension. 



B. Two dimensional treatment 

We now move on to consider the same system as above, but now in strictly two dimensions. 
In 2D the relative coordinate, r = \ri — r^\ y s-state asymptotic scattering wave function is 

V>(r) = = - (cot(5) J (kr) - N (kr)) (r > b) (8) 



given by [9 



where b is the range of the potential and Jo(kr) and No(kr) are, respectively, the zero order 
Bessel and Neumann functions. Before we present our derivation of the 2D pseudopotential, 
it is worthwhile clarifying our convention for the scattering length in two dimensions. 

In ID and 3D systems, the definition of the s-wave scattering length is unambiguous. 
Specifically, in the asymptotic region, and E — > 0, the ID and 3D problems reduce to 



£s»(r) = °> ( 9 ) 



where u(r) = rip(r) in three dimensions and, letting r = \x\, u(r) = ip(r) in one dimension. 
This above equation is solved by u(r) = C(r — a) where a is called the scattering length. 
In other words, a is identified as the intercept of the zero-energy wave function (or its 



extrapolation) on the horizontal r-axis. In 2D, we define u(r) = y/ri[)(r) and obtain the 
asymptotic, reduced radial equation, 

£«M + = ~^{r). do) 

Now, as before, we let k — > and are required to solve the differential equation 

The solution to Eq. (11) is given by the function w(r) = diy/r In r+t^V^ where c?i and d 2 are 
constants of integration. Next, we rewrite our solution in the form u(r) = C\pr (In r — In a). 
Following the same arguments as for the ID and 3D cases, we define the scattering length as 
the node in the asymptotic zero-energy wave function. Note that by definition, the scattering 
length in 2D is strictly positive, which is in stark contrast to the ID and 3D systems where 
a can be of either sign. [l| 

Coming back to the problem at hand, we consider an interaction potential of zero-range 
and utilize the effective range expansion [h]] 

cot(5(k)) = Uln(^) +1 ) , (12) 



7T V V 2 / 

where 7 is the Euler constant, S(k) is the scattering phase shift and a is the previously 
defined scattering length. As before, we can extrapolate our results to a bound state. Thus 
we have i = - (in + 7) which reduces to k = 2ie a . Substitution of the latter expression 
into the scattering wave function, (jSj), gives 

^iH^-M^)) (r>0) - (13) 

The above equation can 



by using the relations 11] 



je re-written in terms of the modified Bessel functions, Iq and K 

J (iy) = l (y) , (14) 

and 

N (iy) = il (y)--K {y) (15) 

7T 

where y is a generic argument. Using Eqs. (fH|) and (TT5T) enables us to obtain the bound 
state wave function, viz., 

tlj(r)=K ( (r>0), (16) 



with binding energy E = — M 4 f \^ . The appropriate 2D pseudo-potential can now be obtained 
by following the same arguments as for the ID case. Namely, we first make the substitution 
r = r(r), and then take the 2D Laplacian of our bound state wave function. Some simple 
algebra gives 

V*) = ^ + r») U (Z^m) + 4r ' ( f 2 '' US 

•• (17) 

where the primes denote derivatives with respect to the specific arguments. We now make 



use of the following useful properties of K , [11] 



and 

,r(r)\ a 9 f f r(r) 



K' 2e" 7 ^ = — — K 2e" 7 ^ . (19) 

V a J 2e~V(r) dr \ \ a J J v ; 

Equations (18) and (19), along with ip(r) = K ( 2e ^ ), allow us to rewrite Eq. (jTTI) as 

V 2 ^W = + r"(r) - ^ W + ^"^ (r) - (20) 

We now make the important observation that the term in the parentheses on the right hand 
side of Eq. (1201) can be rewritten in terms of the 2D Laplacian acting on ln(r(r)), viz., 



+ r »( r ) _ = r (r)V 2 (ln(r(r))) . (21) 

r r(r) 

Use of the expression above, and replacing r(r) by r, allows us to write Eq. (1201 as 

h 2 h 2 d Ah 2 

- F VVM + -V 2 ln(r)r-^(r) = -^^-^(r). (22) 

Finally, substituting £ = and V 2 (ln(r)) = 27r<5 M (f) (see Ref. [8j) into Eq. (HID 

yields 

- ^V 2 ^(r) + ^5 2d (r)r^(r) = E^r). (23) 



Evidently, the pseudopotential 



V™(r) = ^<"(0r|, (24) 



will reproduce the bound state wave function and binding energy for any zero-range potential 
in two dimensions. Notice that in contrast to ID, the Dirac delta function is modified by 



rd/dr, which yields well defined behaviour at the origin. Indeed, this modification is the 
so-called regularization operator referred to in the literature, which is invoked to ensure 
self-adjoint property of the two-body Hamiltonian. (5], 0,0, 12, 3] In fact, our Eq. (24) 
is simply a member of a family of 2D pseudopotentials, [l4j any of which will reproduce 
the shape independent results we have just obtained. In order to clarify this point, let us 
consider the 2D zero-range pseudopotential as derived by Wodkiewicz in Ref. Q, which is 



a particular member of the 2D family: 
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V w (r) 



-a 2 S 2d (f) 



In 



r 2 \ d 



(25) 



where a 2 is a coupling constant, 7 is the Euler constant and L is a characteristic length. 
This pseudopotential is clearly of a different form from that obtained in Eq. (24). We now 
proceed to investigate the effects of Eq. (24) and Eq. (125]) as r — > + . Allowing Eq. (1241) to 
operate on the bound state wave function gives 



V(r)ijj(r) 



27rh 2 



5 2d (r)r 



2e" 7 



2e" 7 r 



As r — > + the first order, modified Bessel function approaches K\ (y) 
sion in Eq. (126]) reduces to 

V(r)^(r) = -^V(r)- 
We now similarly investigate V , and easily obtain 

'2e _7 r\ . / 7 \ / 2e~ 7 



(26) 

^ and our expres- 
(27) 



V w {r)i){r) = -a 2 S 2d {r) 



2(1/ 



K 



-in (y^ 



r 7 
— e 2 ) r 
L 



2e 



(28) 



Given that r 



— 1 for r — > + , we also have the relation Kq 



— In - — - — 7. These two properties allow us to write the small r behaviour of ( 128]) as 



V w (r)i)(r) = -a 2 5 2d (r) In 



L 



(29) 



The coupling constant a 2 can now be related to our scattering length a as follows. In Ref. [3] , 
the bound state energy in terms of a 2 is given by 



E 



4h 2 ir 4ft 2 vr 

-g Afn 2 



ML 2 



(30) 



This energy must be the same as our binding energy given in terms of the scattering length, 



viz., E 



4h 2 



Ma 2 e 2 t ■ 



Equating these expressions gives us the desired relationship between a 2 



and the scattering length a: 



« 2 = ^rr ln ^7— • (31) 



M 

Using this expression along with fl29l) gives us 

U»^(r) = -^ M (f). (32) 

Consequently, we see that as r — > + both our pseudopotential, V(r), and the pseudopo- 
tential of Wodkiewicz, V (r), are operationally equivalent over the appropriate space of 

n 

wave functions. [14J Thus, although the pseudopotentials given by Eqs. (24) and (25) have 
different forms, they both lead to the same low-energy physics. Indeed, in the next section, 
we expand upon this result by showing that the energy spectrum of the ID and 2D systems 
are independent of the details of the interaction potential, provided the range of the potential 
is much smaller than the oscillator length. 



III. ENERGY SPECTRUM OF THE TWO-BODY PROBLEM 
A. One dimensional treatment 

Consider the problem of two identical bosons, each of mass M, confined by a ID harmonic 
oscillator potential. To begin, we will first consider the two particles to be non-interacting. 
Each particle is subject to the potential ^Mw 2 x 2 . In the center of mass, X, and relative 
coordinate, x, we have the Hamiltonian 

H = -— + VV - — + -M c u 2 X 2 , (33) 
2fidx 2 2 P 2M c dX 2 2 ' v ; 

where \i = M/2 is the reduced mass and M c = 2M is the total mass of the system. As we 

will ultimately be interested in the two-body interaction between the atoms, we focus on 

the relative coordinate, where the time-independent Schrodinger equation is given by 



h 2 d 2 iP(x) 
M dx 2 ' A 

2E 



+ -Mw 2 x 2 xl)(x) = Ei/)(x), (34) 



and E is the relative energy. Equation (l34"j) is simplified by defining r] = j— , z = 
I 2 = jj-, which correspond to the dimensionless energy, dimensionless length, and oscillator 
length, respectively. The preceding substitutions transform Eq. ( 1341) into 



d 2 i) 



dz 2 
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+ z 2 i/j = rjip. (35) 



We now assume the solution to fl35l) is of the form ip = e~~zf(z), where f(z) is some 
function of z. By defining y = z 2 and f(z) = w(y) we can write ip = e~zw(y), which upon 
substitution into Eq. fl35l) gives 

d 2 f 1 \ d 1 — r] 



This equation is of the confluent hypergeometric type 



ll|] and as such is solved by a linear 



combination of confluent hypergeometric functions. Equation (36) is of the general form 

d 2 d 

—v{y) + {b-y)—v{y)-av{y) = 0, (37) 



and for non-integral b has the solution [11] 

v(y) = Cl M{a, b, y) + c 2 y x - h M[a - b + 1, 2 - b, y), (38) 



where M is the confluent hypergeometric function of the first kind. H[ The solution to 
Eq. (36) is therefore given by 

™(v) = c X M (l^L, i y) + c^M (^L, ~, y) . (39) 

Equation f[3"§j) . along with our prior definitions gives us the wave function 



^W = ^ 1 M^^,-,^j+c 2 (^)^M^^,-,^jj e ^, (40) 

where we recall z = The ratio, can be extracted by investigating the large z 

behaviour of Eq. (140p . For large argument, y, the behaviour of M is given by 

M(p,q,y)^^y p -^. (41) 

Use of relation (1411) gives us the large 2 behaviour of our solution, Eq. (IIP)) , namely, 



r(i) r( 



3- 



l+n - - 



<'(:)- l c iT^^ + c 2f7^)J 2 "^ e% - ( 42 ) 

For V ; (- 2 ) n °t to diverge at large z, the term in the parentheses must vanish. Forcing this 
term to vanish leaves us with the relation 

* = -ISL (43) 
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ip(x) + -Muo 2 x 2 + V s (x) ip(x) = Eip(x), (44) 



Next, we consider the same system, except now the two particles interact via a short- 
range symmetric potential. We again solely concern ourselves with the relative coordinate 
for which the time-independent Schrodinger equation reads 

h 2 d 2 .. . fi 

~ Mdx~ 2 ^ x) + {I' 

and V s (x) is a generic short-range interaction potential. If we now exclusively consider the 
region where \x\ — > + the harmonic potential vanishes. We are then left with solving the 
problem of a short range interaction potential, which we will take to be of zero-range, so 
that the wave function is given by Equation (3). Note that even if the interaction has a 
finite range, b, the harmonic potential can be ignored if fiu 2 b 2 <C hu, which gives b/l -C 1 for 
the validity of the spectrum derived below. In the the small \x\ region, we have, to within 
a constant 

ip(x) ~ (a - \x\). (45) 

This solution must join smoothly with our solution to the harmonic potential problem. We 
therefore must investigate the small \x\ behaviour of equation Eq. (1401) in an effort to join 
it smoothly with Eq. (145p . For small values of y the confluent hypergeometric function 
M(p, q, y) goes to unity and thus the small z behaviour of Eq. ( 1401) is given by, 



(ci + caCz 2 )*). (46) 
Recalling that (z 2 )^ = we can re- write Eq. (|46j) as 

ip(x) ~ (^-^-V2l- \xfj . (47) 
Relating equations ( )45i) and ( 1471) gives us 



a ci 



which, along with Eq. (43) yields 



a 
1 



(49) 




Equation ( 149]) is identical to the result obtained by Busch et al. [5( but has been derived 
here with no explicit mention of the form of the interaction. Therefore, the energy spectrum 
in ID is universal, and independent of the shape of the short-range interaction potential, 
provided the range is much smaller than the oscillator length. 
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B. Two dimensional treatment 



While the procedure to determine the energy spectrum in 2D closely follows the ID 
treatment, there are additional complications associated with 2D which warrant further 
discussion. Other than dimensionality, the system we consider is identical to that presented 
in section IIIA. Therefore, we immediately write the non-interacting Schrodinger equation 
in the relative coordinate as 

- S^ M(r) + (l MwV - Si^) MW = EM(r) ' (50) 

where u(r) = ^tp(r). Upon making the same substitutions as in the ID problem, we obtain 
the following differential equation: 

-^Lu+( Z *- 7 L)u = r ] u. (51) 

Assuming that the solution is of the form u = \fze~~ f(z) we can rewrite equation floTT) as 

yj 1 w(y) + (l-y)—w(y)-^-+w(y) = 0. (52) 

From a comparison with Eq. (37), we see that this E q. (1521) is also of the confluent hyper- 
geometric kind. The solution to Eq. (1521 . is given by [161] 

w(y) = ciM (^-, l,y\+ c 2 W (^-, 1, yj , (53) 

where M and W are confluent hypergeometric functions of the first and second kind, respec- 
tively. Equation (l5"3l typifies the central complication associated with the 2D case, namely, 
the second argument of W is of integral value; in ID, the second argument of W is 1/2 



When the second argument of W is non-integral, it may be written in terms of M, viz., 16] 



W(a, b, y) = y l - b M{a - b + 1, 2 - b, y), (54) 

which explains why Eq. (39) can be written soley in terms of M. In 2D the integral value 
of the second argument of W does not grant us the ability use Eq. (54). This is not an 
insurmountable obstacle, as we now show. 

Our solution ( 1531 along with the ansatz that u = ^fze~~ f(z), gives 

u{z) = Um f^p, 1, ^ + c 2 W (^-, 1, z 2 ^j ) v^e"^, (55) 
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where we recall z = The ratio, ^, is once again obtained by investigating the large 



behaviour of Eq. fl55|) . where the functions M and W behave as follows, 



16| 



and 



M(p,q,y) -EM,r«e» , (56) 

W(p,q,y) ^ vrcot(Trp)^^-^. (57) 
Use of relations ( 1561) and ( 1571) gives us the large z behaviour of u(z), viz., 

2 

u(z) ~ (ci + c 2 7Tcot(7rp))z 2p ~ 3/2 e^, (58) 

where p = ^j 2 . Since Eq. (jSHD must correspond to a physical state, the term in the paren- 
theses must vanish at large z, leaving 

— = — 7T cot (717?) = ijj(p) — — p), (59) 
c 2 

where ip(p) is the digamma function and the second equality is a fundamental property of 



the digamma function. 



In the presence of a generic central short-range potential, V s (r), Eq. fl50l) reads 

& d 2 / \ /l 2 2 . t r / \ 1 \ 



M rfr^ (r) + + K(r) " M 4^ J M(r) = EM(r) ' (60) 

where again u(r) = y/rip(r). As before, we now consider the region r — > + , in which the 
harmonic potential vanishes leaving us with the problem of a short range potential. For a 
zero range interaction (for finite range, we again require b/l <C 1), V s (r) — > for all r ^ 
and in the asymptotic region Eq. ( 1601) is simply 

d 2 1 

^ w ( r ) + ^ M ( r ) = -fcVr). (61) 

The solution to Eq. (61) has already been given in Eq. (jHJ), which we recall here for conve- 
nience 

u{r) = | (cot(«J) J (kr) - N (kr)) yfi. (62) 

As we are considering the region where r — > + we must explore the small r behaviour of 
(I62p . For y — > + the zero order Bessel and Neumann functions are 

BO 

My) -> 1- (63) 
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and 

N {y)^-{ln{y/2)+j), (64) 

7T 

where 7 is the Euler constant. Use of Eqs. ( |63|) and ( |64l) along with Eq. (12) gives the 
r — > + behaviour of w(r), 

■u(r) ~ (In r — In a) \fr. (65) 

The above solution must join smoothly with our solution for the harmonic oscillator. We 
therefore proceed to investigate the small r behaviour of our solution of the harmonic prob- 
lem, Equation (155|) . For small r, and hence small z, M —* 1. The problem now is with the 
small r behaviour of W. Specifically, what is at issue here is the lack of literature dealing 
with the small r behaviour of W when its second argument is integral. The primary rea- 
son for this void is likely due to the fact that there is no elementary relationship between 
M(p,b,y) and W(p,b,y) when b is non-integral, as in e.g., Eq. (54) Fortunately, in a little 
known paper published more than 70 years ago, W. J. Archibald [17j has developed a useful 
expression for W(p, b, y) for integral second argument, which we present here for b = 1 

W(p, 1, y) = M(p, 1, y)(\n(y) + ftl - p) + 2 7 ) + £) ( f ^^^ [ ) V B , (66) 

where B n = (i + + ... + p+ ^_J - 2(1 + | + ... + £), T(-) is the gamma function and 
again ip(-) is the digamma function. From this expression it is quite evident that the small 
y behaviour will be 

W(p, 1, y) - \n{y) + ^(1 - p) + 2 7 , (67) 
and we find the small z behaviour of ([33]) to be 

u{z) ~ (ci + c 2 (ln(^ 2 ) + ^(1 - p) + 2 7 ))v^. (68) 

Equation (|68|) . upon the replacement of z = -Jj=- can be written, to within a constant, as 

u(r) ~ ^ + ln(r) - hx(l) - \ ln(2) + t^-A + 7 j ^ . (69) 
Utilizing Eqs. (165)) and (l6"9"j) gives 

|- - In(/) - 1 ln(2) + + 7 = - ln(a). (70) 
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Equation (ITOjl simplifies to 



ri + - p) = In ( A ) + 2 In 2 - 2 7 , (71) 



c 2 v ' V 2 « 2 
which, upon recalling that f- = i>{p) — ^(l — p), gives the relation 

^(p) = ln (^) +21n2 -2 7 . (72) 
The spectrum of the system is finally described by recalling that p = — % + |, whence 



^^-i) = b (^) + 21n2 - 2 ^ < 73 > 



Equation (73) differs from the result obtained by Busch et al. by the last two terms 
on the right hand side. This difference can be traced back to the specific form for the 
effective range expansion, viz., Eq. (12), we have used in this paper. If we use Busch's 



expression 



181 ]. cot(8(k)) = (2/n) ln(fca), the last two terms on the right hand side of Eq. (73) 



disappear, in agreement with Eq. (21) in Ref. [5|. Numerically, the difference between our 
Eq. (73) and Busch's Eq. (21) is unimportant, as it would almost certainly not be resolved 
in experiments. The critical point here is that we have shown that the 2D spectrum is 
independent of the details of the short-range potential. In addition, the complications arising 
from dimensionality are not due to the logarithmic singularities in the pseudopotential - as 
suggested in Ref. [5( - but rather from the logarithmic behaviour of W(p,b,y) for integral 
6=1. Indeed, any short-range potential will yield the spectrum above, provided its range is 
much smaller than the oscillator length. It is nevertheless interesting to note that Eq. (73), 
without specifying a form for the potential, has naturally led to the same coupling constant, 
m ,, MW]- 1 which has been used to characterize the strength of the regularized 2D 
zero-range interaction in earlier investigations. [5|] 



IV. CONCLUSIONS 



In this paper, we examined the two-body problem of harmonically trapped ultra-cold 
atoms in one and two dimensions. We have shown that the energy spectra for both ID and 
2D is universal, in that they are independent of the details of the short-range potential, 
provided the range of the potential is much less than the oscillator length. Furthermore, we 
have illustrated that the concept of a zero-range pseudopotential in low-dimensional systems 
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can be easily understood without having to invoke the advanced mathematical language 
of regularization operators. In contrast to more complicated expressions reported in the 
literature, we have shown that our simple 2D zero-range pseudopotential (i.e., Eq. (24), 
without lo gari thmic singularities in the potential) will yield the same low-energy physics 



.an 



13|, ll4j We anticipate the ID and 2D spectra presented here to be verified by current 



experimental methods. 
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